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! As various types of perturbations tend to drive a geostationary 

satellite away from its prescribed position, occasional orbit correc- 
tions have to be carried out by means of a suitable propulsion 
system. In future geostationary missions, low thrust electric propul- 
sion is likely to be applied for station keeping because of 
considerable mass savings. In this paper a station keeping strategy for 

I electric propulsion’ systems is developed. Both the unconstrained case 
■and the case where thrust operation constraints are present are 
considered and tested by computer simulation of a realistic example. 
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Station Keeping of Geostationary Satellites 
by Electric Propulsion 


Summary 

As various types of perturbations tend to drive a geostationary 
satellite away from its prescribed position, occasional orbit 
corrections have to be carried out by means of a suitable 
propulsion system. In future geostationary missions, low thrust 
electric propulsion is likely to be applied for station keeping 
because of considerable mass savings. In this paper a station 
keeping strategy for electric propulsion systems is developed. 

Both the unconstrained case and the case where thrust operation 
constraints are present are considered and tested by computer 
simulation of a realistic example. 
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1. Introduction 


Zz 


The requirement, that a geostationary satellite must always be 
located vertically above a prescribed fixed point of the earth's 
equator, presupposes an equatorial and circular orbit whose rotational 
time coincides with the rotational duration of the earth,. However, 
since an exact geostationary orbit cannot be realized because of 
unavoidable entry errors and the always acting orbit perturbances, 
the limits of allowable deviations from the prescribed position 
play an important role. Depending on the size of the tolerance 
window, the orbit must be corrected more or less frequently. The 
velocity increments needed for this purpose can be generated with 
the aid of the propulsion systems incorporated in the satellite. 

In general the latter are mounted such that thrusts are possible 
in directions tangential and normal to the orbit independent of 
one another. 

Prom the requirement, on the one hand to maintain the limits 
of the tolerance window, and on the other hand to also save fuel, 
one can derive for a prescribed propulsion system the optimum 
correction strategy which specifies the thrust times. Here one 
obtains for the electrical propulsion systems planned to be deployed 
in the future to an increasing degree, considerable differences 
as compared to the today conventional chemical engines because of 
the low thrust level. The advantage of electrical propulsion lies 
in the considerable weight savings for missions of long duration 
[ 1 ], [ 2 ]. 

On the other hand, the low thrust level has a disadvantage in 
that numerous and long duration thrusts are necessary for station 
keeping, which under certain circumstances disturb the remaining 
mission operation. In addition the high energy requirement can 
lead to the fact that the already low thrust level and thrust 
times must be restricted if the limitedly available energy is needed 
at the same time for other mission purposes C3] • 
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If no such constraints need to he taken into account, the 
required thrust times can he calculated analytically. But also 
in those cases, where constraints exist and are known for some time 
in advance, the optimal correction strategy can he found with the /8 

aid of known optimization methods. 

Since the classical orbital elements for the geostationary 
case are in part undefined, it is useful to introduce another set 
of orhit parameters, which are to he designated as "geostationary 
elements" [4]: 


whereby 


D 

h 

1 

P 

q 

1 


= (n-n )/n 

9 g 

- esin(o)+Q) 
= ecos (oo+Q) 

= sin-j sinQ 
= sin-j cosQ 
= M+fl+co— L , 


n 

a,e,i,Q,co,M 

L 


n 


g 


4 


= average motion of the satellite 
= classical orhit elements 
= right ascension of the geostationary 
required point 

= angular velocity of the rotation of 
the earth 

= geocentric gravitation constant 


The elements D, h, 1, p, q, X are also well defined for 
circular and equatorial orbits as dimensionless quantities and 
disappear altogether for undisturbed, exact geostationary orbits. 
They describe the motion of the satellite around the earth and thus 
also the deviation of its apparent position from the required point. 
D is the drift rate normalized with respect to the rotation of the 
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earth; X is the average deviation from the required point while the 
pairs (l,h) and (q,p) essentially represent the eccentricity- and 
inclination vector. 


As a result of the perturbations the geostationary elements 
are functions of time which can be subdivided into secular and 
periodic terms. If one restricts himself to perturbances of the 
first order, one encounters primarily terms with periods of fractions 
of a day, one moon orbit, and one year. Disturbance terms with 
longer periods, e.g. orbiting times of the apsides and the junctions 
of the lunar orbit can be treated during the station keeping problems 
as secular perturbations,, On the other hand the short-periodical 
perturbations with periods below one day must be accepted as 
unavoidable deviations since, with the presently available propulsion 
systems, it is either ineffective or even impossible to control them. 

The disturbances of the geostationary path are caused primarily 
by the three-axis nature of the earth potential, the gravitation 
from sun and moon, and by the solar radiation pressure. The 
analytical description can be simplified [53 because of the special 
properties of the geostationary orbit compared to the very complex 
general perturbation theory. 

Figures 1 to 3 present examples for the time variation of the 
geostationary elements as the result of orbit perturbations. If 
one requires a tolerance window of +0.1° in length and width, then 
- without taking into account the short-periodic effects - the 
north- south limits from an initially exact geostationary orbit are 
violated in this case after 33 days, the east-west limits already 
after about. .20 days. These times are shortened still more if the 
unavailable short-periodic effects are taken into account. 
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2. Station keeping strategy /lO 

In order to prevent the drifting of the satellite outside of 
the tolerance window, the orbit is corrected from time to time. 

This is done with the aid of a propulsion system which consists of 
several engines and which produce the necessary thrusts. 

The specification of the thrust times for all engines taking 
into account the tolerance window and the requirement for minimum 
fuel consumption is the task of the station keeping strategy* As 
already mentioned in the introduction, it is neither possible nor 
desirable here to compensate for all appearing disturbance 
accelerations. Rather it makes much more sense to correct only the 
long-term acting effects from time to time while the short-period 
perturbations are taken into account on an overall basis by 
corresponding reduction of the window limits. Under "long-term" 
we understand here all secular effects and the periodic perturbations 
with periods of more than 30 days. Thus the orbit corrections to 
be made need to be calculated only from the disturbance effects 
which, for instance, are shown in figures 1 to 3 as solid curves. 

The conventional station keeping makes use of a type of 
"wait-and-see strategy" : By means of continuous orbit determination 

one follows, e.g. , the position of points (q,p) in figure 3 long 
enough until it reaches the tolerance limit - reduced by the value 
of the short-period disturbances - and then corrects elements q 
and p by a thrust directed northward or southward in such a way 
that (q,p) is displaced to the opposite-lying edge of the window. 

Similar methods apply for the east-west corrections where the 
conditions are somewhat more complicated only because the four elements 
D, h, 1, X cannot be corrected independently of one another. If the 
perturbations of the eccentricity are only relatively small, it 
suffices to reverse by means of a tangential thrust the algebraic 
sign of the drift rate D. With correct dosing the point (\,D) then 
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describes a parabola which utilizes only a part of the allowable 
tolerance window, while the remainder is taken up by the eccentricity /ll 
effects which were corrected only partially or not at all C6J. 

However, the narrow tolerance windows of future geostationary 
satellites as well as the considerable perturbations caused by 
radiation pressure on the large solar paddles will also require an 
extensive correction of the eccentricity. In these cases at least 
two separate tangential thrusts are necessary for a complete 
east-west correction. 

For the example shown in figures 1 to 3 a north-south correction 
would have to be made, according to this strategy, about every 60 days 
and an east-west correction every 40 days. 

The relatively large north- south thrusts include, as the result 
of unavailable orientation errors under certain circumstances a 
noticeable non-ealculable east-west component, which makes 
premature east^west corrections necessary. Therefore the station 
keeping strategy is usually arranged such that a few days after 
every north-south correction an east-west correction follows 

The conventional "wait-and-see strategy" cannot be used for 
station keeping with electrical engines because the thrust level, 
lower by a factor of about 1000, is not sufficient to achieve the 
necessary inclination changes of 0.2° with a single thrust. This 
could only be done in numerous small steps and would take too much 
time since the perturbation effects are counteracting. If one assumes 
that the TV satellite weighing 1058 kg and planned for 1983 is 
propelled daily for 6 hours with l6.4 mN , one would require about 
170 days for the desired north-south correction. 

When using electrical propulsion systems it makes no sense 
therefore, to wait until the window limits have been violated, 
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but rather to correct the orhit by small amounts in intervals as 
short as possible. For example, one could carry out the north-south 
station keeping in such a way that the point (q,p) is brought back 
[83 daily to the origin (0,0). However, here one must consider the 
small value of the changes of the orbit plane, caused by the disturb- /12 
ances, of only about 0.002° per day, which under certain circumstances 
can be found by orbit determinations just barely or only with 
considerable errors. Thus by the use of such a "one-day strategy" 
one encounters the danger to correct along with increased fuel 
consumption deviations which in reality do not exist, but which 
are only simulated by orbit determination errors. 

In order to avoid this as far as possible, one can, to be sure, 
retain the daily corrections, but specify them each time only in 
accordance with larger time intervals. 

The orbit is determined continuously and by taking into account 
the daily corrections. If a cycle of duration T has passed, the 
duration of the next cycle and the thrust times are specified such 
that the target orbit is reached nominally at the end of the new 
cycle. The corrections to be made for the elements AE^ are 
obtained from the initial values E^° determined from the orbit 
determination and the disturbance effects 6E. during the next cycle. 

(2) AE ± = E^-E?-SE i , 

whereby E? are the elements of the target orbit. 

The thrust times can be calculated from T and from the AEis 
with the aid of formulas which will be derived in the next section. 

Such a long-duration strategy has the following advantages: 

Observations from a longer duration timespan are available 

which make possible a more accurate orbit determination. 
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The deviations., from the required orbit have grown as the result 
of the orbit model- and execution errors after time T to a 
more distinctly recognizable value. 

The time T can he selected appropriately from the requirement 
whereby the magnitude of the various errors and mission- 
engineering reasons can play a role. 

3- Calculation of the thrust times for uncons t rained station keeping 

For a tolerance window of the order of magnitude of + 0.1° all 
orbit corrections to be considered can be considered as differen- 
tially small and those perturbation equations can be used for their 
determination which assume for the geostationary orbit elements the 
following simplified form QV] : 
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1 
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= -^ b T 
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= V (~ b R cosL+2b rpSinL) 
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~ y (b R sinL+2b T cosL) 

~ 2V b N S "'‘ nb 
= 2V b N COsL 


“ - V b R 


Here V = 3* 07^6^7 km/s is the orbit velocity. 

Since only perturbations of 1. order are of interest, we 
introduced on the right-hand side the values of the exact geostationary 
orbit. The quantities b^ , b b^ are the disturbance acceleration 
components in radial, tangential, and normal direction with respect 
to the orbit. They can also be interpreted as the acceleration 
caused by one engine or by a combination of several engines. If one 


9 



assumes these to he constant and considers that the length L increases 
linearly with time, one obtains the orbit correction resulting from 
K thrusts by integration: 


(4) 


(5) 


( 6 ) 


AD ‘ - I j, b Tk T k 


2 * ' 1 

Ah * VS * ( ' b Rk CO£ V 2b Tk Sinb k )Sin 2 rn: k 

k=1 

7 K 1 

A1 = Vk J, (b Rk SlnL k +2b Tk COSL k )Sll '2 nT k 


1 K 1 

A P ‘ Vn , Z , b Nk SlnL k 3in 2 Iu k 
k=1 


1 K 1 
AQ = ^ k =i b Nk cosL k Sln 2 nT k 


A 4 


(7) AX 0 y ^ b Rk T k 

Here b R ^, b^, b^ are the acceleration components, the duration 

and L v the arithmetic mean from the length at the beginning and at 
k +1. 

the end of the k thrust. AX is the correction for the average 
epoch deviation. The correction of the mean deviation at a time Tg 
after the expiration of all thrusts is obtained as the result of 
k drift rate changes as 

(8) 4A = AX 0 -f ^ b TkV L E -L k ) , 
where L_ denotes the length associated with T „ 
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On the left-hand side of equations (4) to (8) one finds the 
desired corrections while the desired maneuver parameters and re- 
appear on the right-hand side, from which one can determine the 
starting times according to the relation 


is) *; 2k - 1 = t„ u tL k -v e G (t ° )] ? K 

fc 2k 

With t we denote the starting time of the cycle, with the 
eastern length of the geostationary required point, and with 6^ 
the sidereal time of the meridian of Greenwich. On the left side 
the odd indices apply for the starting times and the even ones 
for the shut-off times. 

The values for tw^, b^, and tw are fixed hy the propulsion 
system and the mass of the satellite which, with good approximation, 
can be assumed to be constant. 

The fuel requirement is proportional to the total velocity 
increment /l 5 


K 

(10) AV = I b X, 
k=1 k K 


The b^. values are composed of the amount of the acceleration generated 
by the individual engines, which take part in the thrust. Since 

the engines can frequently not be aligned radially, tangentially, 
or normally but rather must be set at a given angle .( figure 5 ). 
only the projections of the thrusts on these principal directions 
are effective so that the value of the acceleration acting on the 
crtdt is ^>| k +b| k +bg k £ b k • 
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From the system of equations (4) to (8) one can derive the 
followings 

K s 3 is a necessary prerequisite for the correction of all 
6 orbit elements. For an unequivocal solution for and 
at K > 3 additional criteria must be used, such as the 
minimization of the fuel consumption. 

One can dispense with radial acceleration components since all 
corrections achievable thereby can be obtained also - and more 
effectively - with the tangential components. 

Since equation (6) contains- only the b^ values and the remaining 
equations only the b^ (and b^) "^ e north-south corrections 
( APt Aq) can be carried out independent of the east-west 
corrections (AD, Ah, A1 , AX), insofar as the propulsion 
system can also generate the north-south or east-west thrusts 
independently of one another. 

For short thrust times the sine^nv can be replaced by the 
arguments so that on the right-hand side only the velocity 
increments AV^ = b^r^. still appear, such as is 

the case for pulse-shaped corrections. 


North-south corrections: 


The corrections of the orbits are determined solely by equations 
(6). Under the assumption of equal normal components of the thrust 
vectors, that is b Nk = b^ = constant, one can show that the following 
solution is optimal from a fuel consumption standpoint: /l 6 


L k = arctan — + 2(k-1)u k=1,...,K 

( 11 ) 

T k = n arcsin ( KTH^T ^ a P 2+a U 2 ') = ti 
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whereby 


- ^ < L-, < ^ for b N > 0 , 

i 

+ 5- ^ Li < ■^•Tt for b. 7 < 0 

2 2 N 

If thrusts of the same value are possible also in opposite 
directions, then b^ k = (-1)^ ^b^. and 

(12) L k = arctan + (k-l)n k=1 , . . . /K 


with 


Equations (11) and (12) state that for a prescribed north-south 
correction by means of K thrusts one must carry out each time at 
a given length position or two oppositely- lying length positions 
of the satellite thrusts of equal duration in its orbit, if the 
fuel consumption is to be minimized at the same time. If the 
point (q,p) is to be transposed in the direction to the origin 
(0,0), that is ( Aq, Ap) = a(q,p), then according to equations (12) 
and (1) the values coincide with the junction of the initial orbit. 

The value of the normal component of the velocity increment 

(I. 3 ) AV n = ^ l b Nkl T k = '"'n N arcsin K]X^T ^P 3 *^ 3 ' 

decreases for an increasing number of thrusts K with the limiting 
value AV^ = 2Vv'Ap 2 +Aq 2 ' for K -*• °°, as one can recognize by the 
example of figure 4 . Thus one should carry out in principle as 
many thrusts as possible and thus as short as possible. However, /17 
with the relatively high thrust level of a chemical engine (~ 10 N) 
one obtains for a 1000 kg satellite mass the optimum limiting value 
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practically already with a single thrust. On the other hand the 
lower thrust level of electrical engines can lead to the fact that 

for K < K crit 


( 14 ) 


Kb 


N 


< Vn\/Ap 2 +Aq 2 = 


no solution exists. This constraint plays a deciding role in the 
design of the station keeping strategy for electrical propulsion 
systems. 

East-west corrections: 

The correction of the longitudinal deviations from the geostationary 
required point is determined by equations (4), (5)» and (8), whereby 
the radial acceleration components are set as = 0, since all 

orbit changes can be carried out more effectively with the 
tangential components b^. At least K = 2 thrusts are necessary 
to satisfy all four conditions. The general solution of the 
transcendental system is not possible in closed form, but one can 
reach the' goal iteratively if one uses for the first step the 
approximation 


( 15 ) 


SlHrllT, 

2 k 


* 2 nX k 


This simplification can be justified by the fact that the 
east-west corrections necessary for station keeping require only 
a fraction of the velocity requirement for the north-south 
corrections. Thus the burning times also become correspondingly 
shorter and mostly lie below 4 hours, whereby the deviation 
between the sine and the argument is about 5$* 


If one first considers the case K - 2, then the equations (4), 
(5)» and (8) can be transformed after the introduction of 


( 16 ) 


b Tk X k 

V k = 


14 


v 


f 


AL — 



as follows: 


Zl8 


(17) 


Vt+V 2 


v i +2v i v 2 cosAL+v| 



2 

3 


AD 


( Ah 2 +A 1 2 ) 


ViSinLi+V25in(L 1 +AL) _ Ah 

v cosL +v 2 cos (L, +AL) - Al 

(Lg~Li ) V-) + (L_,~Li ~AL) v 2 = - Al • 


The solution of the first two equations (17) furnishes the 
values for and Vg as a function of the still undetermined 
parameters AL, whereby L-^ is obtained from the 3* equation. If 
one has found for a specific value AL° the associated values of 
v-^ 0 , Vg 0 , L-^ 0 , then one also has 

(18) V! = v° , v 2 = v 2 / Li = Li+2miti, AL = AL° + 2m 2 Ti 


solutions which, however, generally do not satisfy the k. equation. 
However, this can be obtained by a variation of AL° by iteration. 
Thus one obtains a 2-parameter group of solutions which, however, 
do not necessarily exist for all pairs of numbers m-^, m^- 


The tangential velocity increment is 


( 19 ) AV t 


| ~ 3 | AD ! ^ for v^v 2 

( W^(Ah 2 +Al 2 )- jq AD 2 (l+COSAL)/sin^AL for y v 


> 0 
< 0 


and becomes a minimum for an opposite sign of v^ for AL° = tt. 
Therefore one should be able to select from the multitude of solu- 
tions with oppositely directed thrusts those with the smallest value 
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of ir- AL°. On the other hand, all solutions with thrusts in the 
same direction require the same amount of fuel since the selection 
must be made in accordance with other criteria. 


For moderate accuracy requirements it is not absolutely necessary 
to satisfy all 4 equations (17)* Thus it is customary, for instance, 
for conventional east-west station keeping to satisfy only the first 
3 conditions with AL° = it, that is,, optimally from a fuel requirement, 
while the last requirement is satisfied only approximately by a 
suitable choice of m-^ and nig. Vice-versa one can also proceed in a /l9 
manner where one solves the last 3 equations (17) with A L° = it and 
select from the resulting amount of solutions that solution which 
least violates the drift requirements. 

Such incomplete corrections can be designated as differential 
Hohmann transitions. They have the advantage that they are always 
possible and in addition are optimal from a fuel consumption standpoint. 
The approximation becomes the better, the larger one is able to 
select the values of m^ and mg which, however, is limited essentially 
by Lg or by the time T^, available for orbit correction. 

Complete and at the same time fuel-optimal east-west corrections 
can be achieved by the addition of a third thrust. 


With b RK = 


0 and 


(20) L k = Lo+nyi 


k = 1,2,3 


one obtains from equation ( 5) by multiplication by cos -L° and sine L° 
at first 

(21) AhcosL°-AlsinL° = 0 

and 

(22) tanL° = • 


16 



Furthermore one obtains as an expansion of equation (17) the linear 
system of equations 


\ 

1 1 1 1 


\ 

Vi 


1 A _ \ 

- j iD 

(-I)" 11 (-l)™ 3 

. 

v 2 

- 

y v^Ah^+Al 2 ' 

^L £ -L n L £ -L 2 L £ -L 3 y 

' 

l V3 ; 


1 3 AA j 


which can be solved for v^ for all combinations of m^., for which the 

determinant does not vanish. Thus one obtains a 3-parameter group /2C 

of solutions. The criterion of the smallest velocity increment 

here is not sufficient for an equivocal determination of m^. since 

many fuel-equal solutions result. As an additional criterion one 

can then use, for instance, the maximum longitudinal deviation 

f X ( which appears up to time T-. 
max -i-j 

Naturally the presented method can be expanded for a greater 
number K > 3 of thrusts, whereby the number of the possible combinations 
of the m^ values is also increased. However, the thus resulting 
amount of solutions is limited by the time T™, since per revolution 
of the satellite at most 2 east-west thrusts come into consideration. 


All solutions thus obtained are only approximate because of the 
assumption (15)- The accuracy can be increased by iteration if 
one increases on the right-hand side of equation (17) or ( 23 ) the 


expression 


•7p/Ah 2 +Al 2 


by the term 


K m k 

(24) I (-1) K v 
k=1 


- 2 ■ b Tk 

1 V k 

2 b Tk 
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which for the iteration is set equal to zero and subsequently 

calculated with the values v^. which v/ere obtained from the preceding 
iteration. 

Since the burning times r, are shortened with increasing value 
th K 

of K, the 0 iteration must be the more accurate the larger one 
selects the number of thrusts K. Even the fuel consumption 
simultaneously becoming lower would suggest that one choose K as 
large as possible. However, since the fuel required for the east-west 
correction is already relatively small, there thus results only 
a minor savings which, because of the increased expenditures during 
calculation and the carrying-out of numerous thrusts would pay 
only in special cases. 

4. Optimum station keeping during operational limitations /21 

The station keeping strategies and methods for calculating the 
thrust times, mentioned in the previous sections, are only useful 
if the calculated thrusts can be carried out without limitations. 

However, if one must figure that for mis si on- engineering reasons 
the propulsion system cannot be deployed or only to a reduced extent 
at certain times, the thrust strategy must be determined in another 
way. 


In the following we assume that such limitations are known 
in advance at least for a correction cycle of duration T^ - T q and 
that they are such that within prescribed time intervals no thrusts 
or only reduced thrusts are possible. Then the problem to achieve 
in spite of such limitations a station keeping as accurate as 
possible with minimum fuel consumption becomes an optimization 
problem which can be formulated as follows: 
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Given: 


T o’ T E 

E° 

1 


En 



K 

■*"2k-l ’ ^2k 

b Rk ,b Tk ,b Nk 

b k 

^k-l’^k 


starting- and final time of a correction cycle 
initial values of the geostationary elements 
target values of geostationary elements at the end 
of the cycle 

element changes during the time interval ( T q , T-g ) as the 
result of long duration disturbances 
number of the intended thrusts 

lower and upper limits of the time intervals during 
which one thrust each can occur. (k = 1,2, >■., K) 
radial, tangential, and normal acceleration components 
of k bb thrusts (k = 1,2, „ . . ,K) 

sum of the acceleration values of the engines taking 

4- V] 

part in the k thrust. (k - 1, ..., K) 

estimated values for the starting- and shut-off times 

(k = 1,2 , . . . ,K) 


We seek the following: 


/22 


b 2k l’^2k starting- and shut-off times of the thrust so that the 
secondary conditions are satisfied and that the cost 
function is minimized (k = 1,2,...K). 


The secondary conditions are: 


(25) T 2k-1 " b 2k- 1 " t 2k T 2k k 


19 



The cost function is defined by the following expression: 

6 

(26) 4> = I G (E -E^ ) +G 7 AV ,, 

i.-i 1 ; 

where are suitably chosen masses, the geostationary elements 
at the end of the correction cycle and 


' " * K 

(27) AV = I b u (t„.-t 
k=1 


k vu 2k u 2k-1 ) 


the total velocity increment proportional to the fuel consumption. 


The E- values are defined by 

(28) 


E. = E9+6E.+AE. 
l ill 


(i 1/.../5) 


or 

E« = El+Eip (T e -T 0 ) +6E 6 +AE S 

and by equations (4) - (9) as functions of the ’ t 2k-l , ^2k* The 
long-term disturbances 5E^ may be approximated within a correction 
cycle by linear or 6E^ by quadratic functions. 

Naturally there are many other possibilities to define the 
cost function and the secondary conditions. Thus it would suggest 
itself, for instance, to introduce the exact attainment of the 
target orbit with - E? = 0 as secondary conditions and merely 
to minimize Ay, However, such a method would fail in those cases 
where the target orbit cannot be achieved at all because of 
operational restrictions. 

The formulation chosen here has the following advantages: 
o There always exists a solution. 
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o 


All secondary conditions are linear 


in the variables t 


2 k- 1 ’ 




o The cost function can he differentiated everywhere. 


o The minimization of the cost function • leads to a compromise 

solution between station keeping accuracy and fuel requirement, 
which can be displaced by a suitable choice of the weighting 
into the one or the other direction. 


o By setting the corresponding weight functions equal to zero one 
can limit the method to pure north-south or east-west station 
keeping. 

The above-formulated optimization problem can be solved with the 
aid of the one of the numerous optimization methods found in the 
literature, for which completed computer programs are also available. 

The estimated values ^ 2 k-l° ’ ^ 2 k° are obtained from the 
analytical solution of the unrestricted problem according to the 
method of section 3 * 

5 * Application example /24 

In this section we shall determine the optimum station keeping 
strategy for a randomly selected example. Here we make use of a 
combination of the analytical procedures according to section 3 and 
the optimization methods according to section 4 . First the 
unrestricted problem is solved. From the many solutions we select 
one which is affected the least by the operational constraints. 

It provides the estimated values for the following optimization 
method "Generalized Reduced Gradient Method" [9] , which finally 
provides the desired switching times for the restricted problem. 

The geostationary elements of the initial orbit without 
short-term effects should be on 1/1/1983, Oh WZ: 
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D = -5.664826 • 10 ^ 
h=l = p = q = X = 0. 

The value of D is chosen such that with the secular perturbance of 
X it attains the value of zero. 

The target orbit should be equal to the starting orbit and 
should be reached after 10 days so that the following disturbances 
occurring during the 10 days must be compensated by orbit corrections. 


AD = -11.33 * 
Ah = 18.21 • 

Ai = 59.30 • 


H 

O 

l 

CK 

Ap 

- 6 


10 

Aq 

MO 

1 

O 

H 

A\ 


268.44 • 10~ 6 

- 69.37 * 10 -6 
272.79 1 io “ 6 


These values were calculated according to £ for a satellite 

0 2 

in a 19° west required position with 1058 kg mass and 35 m effective 

cross-section. 

Let the propulsion system, which is planned for the operational /25 
TV-satellite missions, consist of 4 electrical engines with 0.01 N 
thrust each which are attached to the satellites in the sky directions 
NE , NW , SW, SEo Thus one pair of engines must always be started 
simultaneously in order to achieve thrusts in north-south or east- 
west direction (figure 5)- Since thrusts can be produced toward the 
north as well as toward the south, two north-south corrections 
are possibly daily. 

The results shown in figures 6 to 14 were obtained with the aid 
of a computer program which consists of 2 parts. In the 1. part we 
programmed essentially the equations (9) to (24) from which the 
solutions for the unrestricted problem are obtained. Then we 
searched according to the criteria of a minimum fuel consumption, 
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lowest longitudinal deviation, and lowest losses- resulting from 

the prohibited times automatically for the most favorable solution a nd 
modified the constraints accordingly. The satellite data, the 
acceleration vectors produced by the propulsion system, and the 
type of east-west corrections (number of thrusts, complete or 
incomplete corrections) as well as the operational constraints 
are fixed by the input while the target orbit is defined by a 
subprogram which can be exchanged at any time. 

In the second part we call out the subprogram GRGA after 
preparation of the input data, whose results are used to simulate 
the time-related course of the geostationary elements. The 
optimization can be carried out in succession for several weightings. 

The number of the weightings and the weighting factors are input. 

The solution of the unrestricted problem in accordance with 
the method described in section 3 is shown in figure 6. The 20 
thrusts for the north- south corrections last daily 2-times 1.8 hours 
and take place shortly after midnight or during the early afternoon 
( world-time) . 

Of the 3 planned east-west thrusts one with only *«0.3 hours 
duration occurs on the evening of the fifth day and two with about 
equal duration (~0.9 or 0.8 hours) on the morning and evening of 
the seventh day. The selection from the large number of possible /26 
solutions was made in accordance with the criteria of minimum fuel 
consumption requirements and maximum longitudinal deviation as 
small as possible during the course of the 10 days. 

If now this solution is limited by the fact that the engines 
cannot be turned on during certain times, then the optimization method 
defined in the 4. section is being used. 

At first we again select from the large number of solutions for the 
east-west corrections of the unrestricted problems the most favorable 
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one, whereby an additional criterion states that the selected 
solution must be affected as little as possible by the constraints- 
Then the thrust times colliding with the prohibited times are 
shortened appropriately or eliminated entirely. The remaining 
thrust times serve as estimated values for the optimization method 
GRGA for the solution of the restricted problem. 

The weighting factors G^ in equation (26) were chosen such 
that the cost function has a minimum for all elements for deviations 
of about 10 ^ . The prohibited times assumed for this example and 
the resulting thrust times are plotted in figure 7* Because of 
the omission of several north-south thrusts the remaining thrust 
times are correspondingly longer than shown in figure 6. Even the 
east -west thrusts are positioned differently because another solution 
of the unrestricted problem, not affected by the prohibition times, 
was used to form the estimated values. 

In figures 8 to 13 the time course of the geostationary elements 
during the 10-day station keeping cycle is simulated for 4 different 
cases. Perturbances with periods below 1 month are not taken 
into account- 

As expected the station keeping is most accurate for the 
unrestricted case while the optimization leads to satisfactory 
results also for the restricted case- However, if one pursues 
in spite of the constraints the strategy of the unrestricted problem, 
one obtains as the result of the omissions considerable deviations 
from the required orbit, which under certain circumstances can be 
greater than without any station keeping. This is to be feared, /27 

for example, whenever an important east-west thrust is omitted 
(figure 13)* 

The total velocity increment for the 10 days of station keeping 
is increased as the result of the operational constraints only 
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insignificantly by 3*5$ from 2.57 m / s to 2.66 m/s. 

For figure 14- we simulated the motion of the satellites with 
all perturbances and corrections in order to demonstrate the effect 
of position keeping with electrical propulsion engines. Because 
of the close sequence of the orbit correction only a very small 
part of the available tolerance window is used. However, it 
must be pointed out that neither orbit model errors nor execution 
errors nor orbit determination errors were taken into account 
in the simulation. 
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7 . Figures 


/3C 



Figure 1 ; Perturbances of the drift rate D and the longitudinal 

deviation X of a geostationary satellite in 19° west geographical 

longitude, starting on l/l/l983» The mass was assumed to he 1058 kg, 

2 

the effective cross-section 35 m • The parabolic curve shows the long 
duration course with time markings in days, the cross-hatched area 
denotes the variation width of the perturbances with periods under 
one month. The left ordinate scale gives the values of the drift rate 

D defined in equation (1) while on the right the to-be-observed daily 

■ • 

longitudinal change X = Dn + X Q is plotted. Inspite of the favorable 
longitudinal position removed a few degrees from an extreme of the 
earth potential the longitudinal deviation reaches already after 
about 3 weeks the value of 0.1°, which is often demanded as 
tolerance limit for future geostationary satellites. 
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: Perturbances of the eccentricity vector (l,h) of a 

2 

geostationary satellite weighing 1058 kg with 35 m 
effective cross-section in 19 ° western longitude, starting 
on l/l/l983 ■ The curve with the time markings denotes the 
long-term orbit changes while the cross-hatched area 
represents the variation width of the perturbances with 
periods under 1 month. After 50 days the eccentricity 
attains the value of about 3 ' 10 which corresponds 
to a swing with 0 . 03 ^° longitudinal amplitude. 
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2 

satellite weighing 1058 kg with 35 ra effective cross-section at 19° 


western, longitude , starting on l/l/l983« The solid curve shows the 
long-term changes while the dashed curve also contains perturbances 
with periods under 1 month. The l4 day period of a perturbance member 
originating from the moon stands out clearly while all other short- 
period terms are unnoticeably small. Here the tolerance limit of 
0.1° inclination is exceeded after about 33 days. 
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Figure 4: The velocity increment Ay^ necessary for an inclination 

correction of 0.1° as function of the number of thrusts K 
for a north-south acceleration b^ = 10*" km/s. The 
limiting value Av^ lies near 5*366 m/s. The required 
correction cannot be achieved with fewer than 20 thrusts. 






north 



I Sud south 

T 


Figure 5 : The electrical propulsion system, consisting of 4 engines, 

planned for the operational TV-satellite missions. The axes of the 
engines must be inclined toward the north- south direction in order 
to avoid contamination of the solar paddles. Thrusts in all 4 
compass directions can be achieved with the arrangement shown 
by paired engine starts. 
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Figure 6 ; Thrust times for station keeping of a geostationary satellite 
weighing 1058 kg with 35 m effective cross-section at 19 ° western 
longitude from January 1-10, 1983* It was assumed that the propulsion 
system shown in fig. 5 can he operated without restrictions. 

The vertically cross-hatched fields represent the north-south thrusts 
whose daily times are pushed back somewhat from day to day and which 
during one year would move away by a whole day. The combination of 
3 east-west thrusts was selected in accordance with the criteria of 
the minimum velocity increment and the minimum longitudinal deviation 
during the cycle. 
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Figure 7: Thrust times for station keeping of a geostationary satellite 
as in figure 6, but by taking into account prohibited times which 
restrict their operation. They were fixed arbitrarily and are 
plotted as open fields. Since the combination of east-west thrusts 
shown in figure 6 must be voided because of the prohibited times, 
we selected a solution which, although somewhat less favorable with 
respect to longitudinal deviation, was not affected by the constraints. 
Compared to the unconstrained case the fuel consumption is increased 
by about 3 > 5 %- 
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Figures 8 - 13 : 


Simulation of the time course of the 6 geostationary /37 
elements during the 10 day station-keeping cycle 
without taking into account perturbance effects 
with periods under one month. 


curves Is 
(curves 2: 
curves 3 J 

curves 4: 


shape without station keeping 
shape for unconstrained station keeping 
shape for operational constraints not 
included in calculations 

shape for included operational constraints 
and optimized station-keeping strategy. 
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Figure 8 : 
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Figure 12: 







Figure 14 ; The area occupied in a tolerance window of +, 0.1° at 19° 
west as the result of longitude and latitude satellite movements 
during the 10 day cycle, taking all perturbances into account. 

dark region: for unrestricted station keeping 

cross-hatched area: for optimized station keeping with constraints 

empty area: without station keeping. 




